Ashkin- Teller criticality and pseudo first-order behavior 
in a frustrated Ising model on the square lattice 
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We study the challenging thermal phase transition to stripe order in the frustrated square-lattice 
Ising model with couplings Ji < (nearest-neighbor, ferromagnetic) and J2 > (second-neighbor, 
antiferromagnetic) for g = J2/IJ1I > 1/2. Using Monte Carlo simulations and known analytical 
results, we demonstrate Ashkin- Teller criticality for g > g* , i.e., the critical exponents vary contin- 
uously between those of the 4-state Potts model a,t g — g* and the Ising model for g — > 00. Thus, 
stripe transitions offer a route to realizing a related class of conformal field theories with conformal 
charge c = 1 and varying exponents. The transition is first-order for g < g* = 0.67 ± 0.01, much 
lower than previously believed, and exhibits pseudo first-order behavior for g* < 5 ^ 1- 
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Breaking a four-fold symmetry (Z4) at a thermal phase 
transition in a two-dimensional (2D) system can lead to 
a very special kind of critical behavior, where the critical 
exponents depend on microscopic details, in contrast to 
the universal exponents at normal transitions. Contin- 
uously varying exponents are possible in conformal field 
theory (CFT) when the conformal charge c > 1 [H, 
A well-known microscopic Z4-symmetric model realizes 
this behavior for c = 1; the Ashkin- Teller (AT) model 
[3, 01 consisting of two square-lattice Ising models con- 
pled by a four-spin interaction. A Z4-symmetric order pa- 
rameter also obtains in the standard ferromagnetic Ising 
model when sufficiently strong frustrating (antiferromag- 
netic) second-neighbor couplings are added, in which case 
the low-temperature ordered state exhibits stripes of al- 
ternating magnetization (which can be arranged in four 
different ways on the square lattice). While c = 1 AT 
criticality has been suspected at the stripe transition, 
it has not been possible to demonstrate this convinc- 
ingly for a wide range of couplings (only at very strong 
second-neighbor coupling, where the system reduces to 
two weakly coupled Ising models 0, ) . It has even been 
difficult to determine whether the transition is contin- 
uous or first-order A convincing demonstration 
of robust c = 1 criticality with varying exponents would 
be important, since stripe ordering can be achieved in 
many different systems [isj and represents a more re- 
alistic route to c = 1 and varying exponents than the 
somewhat artificial AT model. 

Using Monte Carlo (MC) simulations of the frustrated 
Ising and AT models, along with known analytical results 
for the AT model, we here demonstrate c = 1 scaling 
at the stripe transition, with continuously varying expo- 
nents matching those of the AT model. We also uncover 
a remarkable pseudo first-order behavior associated with 
this type of continuous phase transition, which had not 
been previously noticed and which led to wrong conclu- 
sions regarding the nature of the stripe transition. 

The 2D frustrated Ising model with couplings Ji < 



(ferromagnetic) and J2 > (antiferromagnetic) is defined 
by the Hamiltonian 



(1) 



where first and second (diagonal) neighbors on the square 
lattice are denoted by {ij) and {{ij)), respectively, and 
(T; = ±1. A stripe state obtains at low temperature (T) 
when the coupling ratio g = J2/\Ji\ > 1/2. Even though 
this model represents one of the most natural and sim- 
ple extensions of the standard 2D Ising model, its stripe 
transition remains highly controversial. The most basic 
question under debate is whether the transition is contin- 
uous for all g > 1/2, or there is a line of first-order tran- 
sitions up to a point g* . The next issue is to establish the 
universality class of the continuous transitions. Early nu- 
merical and analytic approaches supported the idea that 
the transition is always continuous, but with the criti- 
cal exponents changing with g. Some variational stud- 
ies 14 1 and recent MC studies 0, Q have, however, 
found a line of first-order transitions for 1/2 < g < 1. A 
very recent MC study Q used a double-peak structure in 
energy histograms to conclude that the transition is first- 
order at least up to g = 0.9. For large g the transition 
should be AT-like, based on an effective continuum the- 
ory of two weakly coupled Ising models 0]. For smaller g, 
the AT behavior could either continue down to the point 
at which the transition turns first-order (with the criti- 
cal end-point of the AT line corresponding to the 4-state 
Potts model); or some other behavior that is beyond the 
AT description may apply instead. The MC results of 
Ref. could not distinguish between these scenarios. 

We show here that the stripe transition is first order in 
a much smaller range of couplings than claimed in Ref. 0; 
for 1/2 < g < g* with g* « 0.67. For g > g* it is continu- 
ous and in the AT class. The exponents change continu- 
ously with g as in the AT model \M, with g* correspond- 
ing to the 4-state Potts model and g 00 to stan- 
dard Ising universality [3| • First-order indicators (neces- 
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sary but not sufRcient), e.g., multiple peaks in energy and 
order-parameter distributions, lead to over-estimation of 
the region of discontinuous transitions. As we will show, 
this is because the 4-state Potts model and neighboring 
transitions in the AT model exhibit pseudo first-order 
behavior, though these transitions are known to be con- 
tinuous The Potts point had not been identified in 
previous works, and the pseudo first-order behavior was 
interpreted as actual first-order transitions. 

Stripe transition. — The striped phase is characterized 
by a two-component order parameter (rnx,my) with 
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where {xi , yi ) are the coordinates of site z on an L x L 
periodic lattice and N = L^. We define = + niy 
and the stripe susceptibility x = N{{m'^) — (|77i|)^)/r. 
We employ the standard single-spin Metropolis MC al- 
gorithm and use | Ji| = 1 as the unit of T. 

We analyze the peak value C^a^iL) of the specific heat 
and Xmax(i) of the stripe susceptibility. By standard 
finite-size scaling arguments [l^, Cmax(i) ^ L°'^'^ and 
Xmax(i) L^/" . For first-order 2D transitions these 
quantities should instead diverge as L^. Examples of the 
scaling behavior are shown in Fig [TJ For the L < 256 
systems studied, the exponent a/v, estimated from the 
slope on the log-log scale, decreases with increasing g (it 
is close to 2 for g = 0.5f ), while 7/1/ remains close to 7/4 
for the g- values displayed here (also approaching 2 closer 
to g = 1/2). These behaviors can be affected by scaling 
corrections, and the slopes for moderate sizes may not 
reflect the true exponents. Fig. [T] also shows results for 
the 4-state Potts model, for which it is rigorously known 
that a/v = 1 and ^/v = 7/4, but there are multiplicative 
logarithmic scaling corrections that affect the behavior 
strongly for lattices accessible in MC simulations [l7| . 
For g in the range 0.66 — 0.70, the J1-J2 scaling agrees 
well (up to factors) with that of the Potts model. 

According to CFT [l|, [2[, the exponents can change 
continuously when the ordered phase breaks Z4 symme- 
try. There are two known microscopic scenarios, exempli- 
fied by: (i) The XY model in a four-fold anisotropy field 
114. For /14 = there is a Kosterlitz-Thoulcss transition 
[19l l20ll . while \hi\ -^00 gives standard Ising universal- 
ity [21| . (ii) The line of fixed points connecting the Ising 
and 4-state Potts points in the square-lattice AT model 
Assuming that one of these scenarios applies to the 
continuous transitions in the J1-J2 model, ^/v = 7/4 is 
fixed because this holds always for both (i) and (ii) . This 
explains why x/L^/"* is almost constant for all cases in 
Fig. [ijb). However, the observed variation oi a/v with g 
in Fig.[TJa) (and at larger g) rules out scenario (i), since 
a = for all the fixed points there. Then the natural 
scenario left to consider is that there is a line of fixed 
points corresponding the AT model. 
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FIG. 1: (Color online) Peak value of the specific heat and 
stripe susceptibility vs L for the J1-J2 and 4-state Potts mod- 
els. Factors corresponding the asymptotic Potts scaling have 
been divided out. In (a) the curves through the g = 0.68 and 
Potts data are fits to the Potts form C = aL\n{L/h)~^^'^ . A 
log-correction should also be present in (b) but the expected 
Potts form X = aL^^* ln(L/b)~^/* [13 is not seen for L < 256. 



The AT Hamiltonian can be written as 



H 



TiTj + KaiajTiTj), 



(3) 



where two Ising variables, (Ti, t;, reside on each site i and 
are coupled to each other through K . The ferromagnetic 
phase of the AT model ((ctt) 7^ and (a) = ±(t)) breaks 
Z4 symmetry and the order parameter can again be ex- 
pressed as a vector {m^, m^) with 
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The relevant Z4 transitions take place when K e [0, 1], 
where X = corresponds to two decoupled Ising models 
and K = \ can be mapped to the 4-state Potts model. 
The r > transitions for K e [0, 1] arc all continuous, 
with the exponents depending on K (but = 1/4 always) 

i- 

Binder cumulant. — To further understand the nature 
of the phase transitions, we probe the Binder cumulant 
of the stripe order parameter. For a 2-component vector 
order-parameter, the cumulant is defined as 



[7 = 2 1 



1 (m^) 

2 (^^ 



(5) 



where the factors are chosen to make C/ in the dis- 
ordered phase and U ^ 1 in the ordered phase in the 
thermodynamic limit. We define U for the AT model in 
the same way, by using = + . 

For continuous transitions, the cumulant typically 
grows monotonically upon lowering T and stays bounded 
within [0,1]. It approaches a step function at Tc as 
L — > 00 |22l |. For a first-order transition, it instead 
shows a non-monotonic behavior with T for large sys- 
tems [23[, developing a negative peak which approaches 
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FIG. 2: (Color online) Binder cumulant vs temperature for 
(a,b) the J1-J2 model at g = 0.55 and 0.70, and (c,d) the AT 
model a.t K = 1 (the 4-state Potts point) and K = 0.95. 



Tc and grows narrower and diverges as (in two dimen- 
sions) when L — > (X). The nonmonotonic behavior can be 
traced to the emergence of multiple peaks in the order- 
parameter distribution (reflecting phase coexistence), as 
discussed in the context of the J1-J2 model in (23 |. 

As is clear from Figs. [2ja,b), U{T) of the J1-J2 model 
indeed develops a negative peak that grows with increas- 
ing L for g = 0.55 and g = 0.70. As g increases, the 
system sizes needed to observe a peak also increase, indi- 
cating a weakening discontinuity of m. The dependence 
of the peak value C/min on L is shown in Fig. [3] for sev- 
eral values of g (where C/min > corresponds to a local 
minimum). Interpolating these data, we can extract a 
length Lo{g) where C/min crosses zero for given g. Lo{g) 
grows with increasing g and diverges when g w 0.82 (i.e., 
for larger g there is no negative peak). The vanishing 
of the negative peak might be taken as an estimate of 
the location g* of the multicritical point. One could 
also examine the L{g) at which a minimum in C/(T) first 
forms (but is not yet negative). This gives a still higher 
value of g*. However, such procedures, or ones based on 
multiple peaks in the order-parameter or energy distri- 
bution, over-estimate g* because these features can ap- 
pear also for continuous transitions. Indeed, as shown in 
Figs [21Ic,d) and [3l we observe negative cumulant peaks 
also in the Potts and AT models, in spite of these models 
being rigorously known to have continuous transitions. 

The most natural scenario suggested by these data is 
again that the J1-J2 model at a point g = g* corresponds 
to the 4-state Potts model. It therefore exhibits pseudo 
first-order behavior for some range of g- values above g*, 
just like the AT model does for K at and slightly be- 
low the Potts point K = 1. The very similar behavior 
of Cmax(i) for g = 0.68 and the 4-state Potts model 
in Fig. [IJa) already suggests g* w 0.68. Below we will 
present further evidence of g* = 0.67 ±0.01 indeed being 
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FIG. 3: (Color online) Peak value of the Binder cumulant vs L 
for the J1-J2 model at different g values and the 4-state Potts 
model. For the latter (1) in the legend refers to the original 
data and (2) has rescaled size, L — >■ aL, with a ~ 2.36. 

a 4-state Potts-universal point. 

We presume that the stripe transition is first-order for 
g < g* , instead of some alternative (and unlikely) ex- 
otic behavior outside known scenarios for Z4 symmetry- 
breaking. The discontinuities are always very weak, how- 
ever, so that the asymptotic first-order scaling behavior 
cannot be observed in practice. For instance, the nega- 
tive U{T) peak for g < g* and the specific heat appear 
to diverge much slower than the expected form. Sim- 
ilar behavior can also be observed for the 5-state Potts 
model, which is a well-known prototypical example of 
a weak first-order transition [l6j with anomalously large 
correlation length ^ [2^ . To observe clear first-order scal- 
ing, one has to use lattices with L » ^. 

Histograms — We next examine in some more detail the 
pseudo first-order signals in the AT and J1-J2 models, 
using the probability distribution of the order parame- 
ter and the energy. Some aspects of the full distribu- 
tion P{mx, my) were discussed in 2^. Here we consider 
P{irL^). It is well known that phase coexistence at a 
first-order transition leads to a double-peak distribution 
in a narrow window (of size in two dimensions) 

around Tc- For L 00, the distribution approaches two 
delta-functions (at to = and the value m of the order- 
parameter just below Tc), with weight shifting between 
the two across the narrow Tc window. A double peak in 
the energy distribution corresponds to latent heat. 

As would be expected based on the negative Binder cu- 
mulants, the pseudo first-order behavior also gives rise to 
double peaks in the order-parameter and energy distribu- 
tions close to the transition. This is shown in Figs.|4|a,b) 
for the 4-state Potts model. Similar behavior can be 
seen in the J1-J2 model, for which results are shown in 
Figs.[4Uc,d) at g = 0.67 « g* (where the transition should 
be weakly first-order) . In Ref. Q a double-peaked energy 
histogram was seen all the way up to 5 = 0.9 on large 
lattices. This was taken as a first-order transition, while 
our results show that this is merely pseudo first-order 
behavior [2^. Combining the results, we conclude that 
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FIG. 4: (Color online) Histograms of the order parameter 
and the energy E for (a),(b) the 4-state Potts model (the 
= 1 AT model) and (c),(d) for the J1-J2 model at g = 0.67. 
Here T is very close to Tc, chosen such that the two peaks in 
the energy histograms are of the same height; for the Potts 
model T/K = 3.64231 for L = 128 and 3.64460 for L = 256, 
while for the J1-J2 model T/Ji = 1.2014 for L = 128 and 
1.2004 for L = 256. 



there is pseudo first-order behavior in the J1-J2 model 
from the Potts point g* « 0.67 up to at least g = 0.9. 

Potts point. — If the continuous transitions indeed be- 
long to the AT class for g G [5*, 00], then a way to es- 
timate the Potts point g* more precisely is to use the 
universal Binder crossing value U* (sec Fig. [2]) of curves 
U{T) for different L at fixed g. For the 4-statc Potts 
model (the K = 1 AT model), we estimate U* = 0.792(4) 
by extracting crossing points between data for pairs 
(L, 2L) and extrapolating to L = 00. Examples of the 
finite-size scaling of crossing points are shown in Fig[5ja). 
In the J1-J2 model, U*{g) increases monotonically with 
(7, as shown in FiglSjb). We can now estimate g* by 
equating U*{g) to the 4-state Potts value. This gives 
g* = 0.67 ± 0.01, in good agreement with the Potts-like 
behavior of the specific heat shown in Fig. [Ija) and the 
histograms in Fig. 

As a further test, we use a method inspired by the flow- 
gram technique 27[: Since the order-parameter distribu- 
tion should be universal (up to scale factors) at Tc for a 
continuous transition, and the negative peak in U{T) rep- 
resents one aspect of this distribution {U* being another 
one) the value UminiL) should either approach a universal 
value (as U* does) or diverge in a universal manner when 
L ^ 00. Then, UminiL) for the 4-state Potts model and 
the J1-J2 model at g* should, if the models are controlled 
by the same fixed point, collapse onto the same curve for 
large systems, once a rescaling L — > aL is introduced for 
one of the models. This is indeed the case if, and only if, 
g is close to the g* point estimated above. Scaled Potts 
data match best the g = 0.66 curve in Fig. |31 It is not 
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FIG. 5: (Color online ) (a) Binder cumulant crossing points 
for {L,2L) system pairs. The approach to U*{L = 00) is 
expected to be governed by a nonuniversal scaling correction. 
The curves show fits of the form U — a + b/L'^. (b) L — ^ cxd 
extrapolated U* of the J1-J2 model compared with the Potts 
result (the lines indicate the estimate with error bar). 

clear whether t/min here diverges very slowly or converges 
to a finite value. 

Conclusions. — By combining many mutually consis- 
tent signals, we have demonstrated a point g ~ g* ~ 

0. 67 ± 0.01 at which the J1-J2 model is controlled by 
the 4-state Potts fixed point. The scenario of the criti- 
cal curve for g G [g* , 00) being in one-to-one correspon- 
dence with that of the AT model for K £ [1,0) is the 
only possibility consistent with known CFT scenarios and 
MC results. The pseudo first-order behavior uncovered 
here, along with the very weak first-order transitions for 
g € {1/2, g*), is what made it so difficult to correctly 
characterize the nature of the transitions until now. 

With the nature of the transition now understood, it 
will be interestin g to study other aspect of the model, 

1. e., the kinetics j28l [29j at the weakly first-order and 
pseudo first-order transitions. We also note that many 
interesting quantum problems involve the same kind of 
order-parameter symmetry in two dimensions, e.g., stripe 
states in models of high- Tc superconductors 3^, 31 1 and 
valence-bond-solid states in quantum antiferromagnets 
[3^ . The T > transitions in these systems, as weh, 
may be impacted by the issues we have pointed out here. 
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